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1. INTRODUCTION 


This report documents the results of a six-month study on the pointing 
control of a microwave antenna of the Satellite Power System (SPS). The study, 
performed under NASA contract NAS8-33604, may be considered as the continuation 
of an early effort on the accuracy analysis of pointing control system of 
SPS, which was conducted under NASA contract NAS8-33065. Primary goals of 
the study encompass three areas, namely, the SPS antenna pointing error sensing 
method, a rigid body pointing control design study, and approaches for modeling 
the flexible body characteristics of the solar collector. 

Accuracy requirements for the antenna pointing control consist of a 
mechanical pointing control accuracy of three (3) arc-minutes and an electronic 
phased array pointing accuracy of three (3) arc-seconds. Results of the 
study, based on the factors considered in current analysis, show that the 
three arc-minute overall pointing control accuracy can be achieved in practice. 

Section 2 of this report deals with the antenna pointing error sensing 
method and its performance. Among the topics discussed are the pointing error 
sensing method, received signals, subarray signal processing, array signal 
processing, effect of noise, and effect of vibrations. 

Section 3 presents the details of a mechanical pointing control design 
for the antenna. Included in the section are the baseline satellite con- 
figuration and a set of baseline parameters, an estimation of motor rating 
and motor parameters, slip-ring friction characteristics, the effect of 
centrifugal force-induced torque on antenna due to gimbal rotation, the 
effect of variation of antenna moment of inertia due to its changing attitude, 
control system modeling, control system design, system simulation, and the 
effect of noise in sensed pointing error signal. Also included in the section 
is the effect of collector vibration on the accuracy of pointing control. 
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Shortly after the start of the study, it was found that the required 
effort for an analytical modeling of the flexible structure of the solar 
collector was immense, way beyond the planned effort. As a result, work 
in this area was redirected into two directions. The first direction is 
to explore the approaches for obtaining an analytic flexible model, and 
the second direction is to estimate the effect of the collector's flexible 
body vibration on the antenna pointing control accuracy. Section 3 is 
devoted to the discussion of the approaches for the flexible collector 
modeling . 

A list of references and two appendices are included for the con- 


venience of readers. 
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2. ERROR SENSING METHOD FOR ANTENNA POINTING 

In the final report of an earlier effort [1] a method was proposed 
for developing angle error signals to be used by angle servomechanisms 
to point the antenna toward the ground station. In this effort, the 
proposed method has been shown to be feasible. In addition, a normalization 
technique, similar to that used in angle tracking radar, has been added 
that produces error signals that are independent of received signal magnitude. 
Finally, some results related to affects of vibration and noise have been 
derived. 

2 . 1 Sensing Method 

The basic concept is illustrated in Figure 1. Two intersecting linear 
arrays formed from subarrays of the main array are used to develop error 
signals. The array normal to axis 1 is used to develop an error signal 
proportional to the error that the antenna pointing direction makes relative 
to the desired direction (tov/ard the earth site) in a plane for which axis 

1 is normal (error about axis 1). The other array generates a signal 

proportional to the (space) orthogonal angle in the plane for which axis 

2 is normal. The two error signals are used by angle servomechanisms to 

turn the antenna about axes 1 and 2 to cause the errors to approach zero. 

Each array is presumed to have 2N + 1 = M subarrays with the cormon 

subarray at the antenna's center serving as the "reference" for the two 
arrays. In general, the two arrays do not have to have equal numbers of 
subarrays; however, since accuracy requirements are likely to be the same 
in both coordinates, there seems to be no reason to not make the two 
identical, as is assumed in following discussions. 




to coordinate 2 


Subarrays applicable 
to coordinate 1 


Figure 1. Subarrays applicable to gimbal servo error generation. 



2.2 Received Signals 
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We define antenna-centered coordinates as shown in Figure 2 where 
X and y are axes in the plane of the antenna's face. Axis z is the 
antenna's true pointing direction. A wave arriving from the earth site 
makes an angle e from the z (boresight axis) direction; it has an 
"orientation" defined by an angle <p in the xy plane. Using these definitions 
we find the signal received by a typical subarray located along the x 
direc.tion. Only the x axis array is analyzed since the other behaves in 
a simi lar manner. 

The error angles sensed by the system in Figure 2 are 0 (in the xz 

X 

plane) and (in the yz plane). The various angles are related by 

tan(e^) = tan(e) cos(cj)) (1 ) 

tan(e^) = tan(e) sin((^>) (2) 

from geometry. 

A typical subarray pattern in the SPS system will have a uniform 
illumination function, have a rectangular shape (side-lengths X and Y) and 
a one-way voltage pattern, denoted G(e, 4 )) by 

G(e,({)) = Sa sin(e )cos((^))]Sa[^ sin( 0 )sin((j))] (3) 

A A 

where 

Sa(c) = (4) 

and A is wavelenght. If we now recognize that |0|<<! for any reasonable 
SPS system, then |9j^l<<l and |6y|<<l, so (1) through (3) give 

0j^ ^ 0 cos((j)) (5) 


0 ^ 0 sin ^ 


( 6 ) 



6 





Figure 2. Antenna-centered 


coordi nates . 
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G(6,4>) ^ Sa[(TrX/x )0 cos((j>)] Sa[(7rY/x)9 sin((j))] 

= SaCirXa /A]Sa[iTY9 /a] = G(e ,9 ) (7) 

X y ^ y 

\ 

Now consider the array of subarrays along the x axis as shown in 
Figure 3. Subarray locations along x are nX, -N<n<N, and they have dis- 
placement in the z direction of and tilt angles (tilt in the other 
coordinate is denotes the voltage produced by the nth sub- 

array by a wave having error angles e and e , then 

X y 

o -w 2irz 

V = Ap, G(0 +0 , 9 +9 )cos[(|> .+ sin(e )+ — ^ cos(9 )] (8) 

n 0 ^ X xn y yn ^ ‘-^pilot A ^ x A ^ x^*" ' ' 

where Ag is the voltage produced when the subarray points toward the 
pilot source without tilt error and ^p-j-jQ^ is the phase of the arriving 
pilot signal at the reference (n=0) subarray. Since 10^1<<1, 

2irnX9 2ttZ 

Vn % Ag G(e^+9^^, ^ + ^]. (9) 

2.3 Subarray Signal Processing 

We next show that desirable error sensing occurs if the signal pro- 
cessing of Figure 4(a) is available at each subarray in the line array. 
Straightforward analysis gives 


(+)] 


(-) 


, A 2TTnXe 27r(z -z^) 


27TnX0 2ir{z-Zf.) 


(10) 
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Figure 3. Subarray locations. 


9 








10 


The signals e and e are assumed conveyed to a central processor for 
^ n n 

the array, as shown in Figure 4(b). 

2.4 Array Signal Processing - No Errors 

Suppose no tilt or displacement errors are present. Then 

9 =0, and z =0 for all n. Now (10) becomes 

yn ’ n 


(+) 


(-) 


The signals e 


(+) 



and e 


(-) 


G(9^, 9y) V®q 

+cos[n(M 9^-a^ 

in Figure 4(b) become 


)] 

)]}. 


(H) 


,(+) 

,(-) 


) = f G(9^. ep ( 


sin[j(^ 8^+6q)] 


where 


- sin[^(M 9^.9 )] 


M=2N+1 


(12) 


(13) 


and we have used the known series 


N 


sin[^4^ B] 


I cos(ns) = 3 

n=-N sin(“/ 2 ) 


( 14 ) 



Next, define a subarray beamwidth (in rad) between -392dB points 


n 


as given by 


and 


sub 


= x/x. 



9 u 
sub 


Sq/2.. 


Then (12 becomes 


(15) 

(16) 


(+) 


1~ ( 9 , 0 ) { 

_lJ 2 ' x’ y' 


(-) 


sin[^(e +0- )] 
^ub " 1 

s1n[^ (e +e J] 


sub 


X q 


sin[^l!— (e -9„)] 


sub 


X q 


siiCl (0v-0o)] 


sub 


X q 


(17) 


Now if we define an "array pattern" G (8 ) by 

a X 


sin[ 


Mtt 


G (e ) = 
a ^ X ' 


®sub ^ 


M 7T 9 

9^] sin(:; -) 


■9 


TT 9 

M sin (- — -) 
®sub 




sub 


M IT 9 


(- 


for [9 |<<9 ,/tt 
‘ x' sub' 


9 


sub 


(18) 


then 


e^"^' A M 

^_^}=^S(0x, 8y) (G3(9x+Bq) 


( 19 ) 


A sketch of G^(s^) is shown in Figure 5 for M=3,5, and 9. The function 
is symmetric about z=0 and + tt and about z=+tt/2. Thus, there is another 
peak in the function at z=+tt. For M large G (a ) approaches a sin(x)/x 
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function for small magnitudes of x. The responses (peaks) at +tt never 
show up in the response (17), however, because G(e , e„)=0 for z=ttX8 /x = 

A jr A 

TT0 /e . = 7T from (7). 

X sub ' 

Next, define array sum and difference patterns by 

Gz(Sx) - 1 ^ (20) 

respectively. Equations of (19) become 

( 22 ) 

(23) 

The response e^ in Figure 4(b) becomes 

(24) 

The behaviors of the sum, difference, and "normalized difference 
pattern" S^((^^)/Gj;(ej^) are shown in Figures 6-9. In these figures we 
define an "array basic beamwidth" by 

®B ■ ®sub''''' ° 

Figure 6 shows that the sum pattern can be of the usual single-peak form 
or double-peaked. The transition takes place when 9q/0g = 0.6625. Figure 
7 shows typical difference channel responses with Figure 8 showing that the 
slope of the response at the origin is maximum when 9q/9g about 0.66. 
Figure 8 also shows the loss in gain of the sum channel as a function of 







Figure 6. Sum pattern. 






Figure 9 . Normalized difference pattern 
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Fq/eg. Finally, Figure 9 shows the normalized difference pattern. From 
these data, which apply for M larger than about 9, it is clear that 
near 0.5 should be a good system choice based on small sum channel gain 
loss and beamwidth broadening, large error slope, and linearity of the 
normalized difference pattern. 

For ? /e„ =0.5 the slope of the difference pattern on boresight is 
q b 


d G fe ) 

T A X 

slope = 


-4 




-4 M 


^ ^ ®sub 


(26) 


The slope of the normalized difference pattern on boresight is 


slope 


d 9. 


e,(e ) 

r_9_iL.T 


9 =0 
X 



-2 M 
® sub 


(27) 


2 . 5 Effects of Noise 

Assume each subarray output contains a bandpass gaussian noise N^(t). 

Such noises can be described by [2] 

N,(t) = t] - s1n[.p.,^P t] (28) 

where '"^gpj(t) and independent, gaussian random variables for any 

"t and a given value of n. They are also independent between subarrays 

p p 2 A 2 

(different values of n). Furthermore E[N^ (t)] = (t)] = (t)] = 

is the noise power (assumed the same at each subarray). 

S traightforward analysis of Figure 4(a), assuming the LPF's pass 

the full noise bandwidth, gives the noises present in e^"*"^ and e^ \ 

denoted N and N respectively, 

n n 


(t) = [Nap(t) cos " \n (t) sin(9^)] cos(nep) 


(29) 
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sin(e^) - cos(e^)] sinCne^) (30) 


Here ^l^p-j “ ^pilot ^ ®T’ arbitrary pilot source 

phase angle. 

After summing the various noises as implied in Figure 4(b) the 
noises on the signals e^“^ and e^*^\ denoted N^'^(t) and M'^^(t), respecti vely, 
are found. The powers in these noises follow taking the mean-squared 
values: 


2 M . 1 »a) 


2 ^ rr r.i( + ) 2t _ rf^ . I 'q' i 

Oj. = E[{ 1 'J (t)} ] = [j + 2 sin(8q) ^ 

cr^=E[{N (t)} ] - [^ - 2 sin(e ) ^ ' 


(31) 

(32) 


Figure 10 plots (31) and (32) versus Q^/Qn “ ^ Clearly, for squint 

2 2 

angles of interest (e^ = 0g/2) we have so output noise powers 

are equal : 


a 


I 


2 



M c^//2. 


6q/0B " 


(33) 


2.6 Effect of Vibration (Displacement) Errors 

We assume that there is symmetry in the vibration of the overall 
array so that the geometry of Figure 11 applies. Since subarray positions 
are 

= n X, -N < n < N, (34) 

then 

= z(x^) = Az cos[(2K-l) 7T X n/D]. (35) 

For small displacements Az << D, as is the case in practice, it is readily 
shown that the vibration induced tilt is 

, -az(2K-1)it 
xn D 


sin[(2K-l) 7T X n/D]. 


(36) 
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7 r x/Pj 

D = A Af re. A/A/A V/amet^r 
K ' VtBRATtoN Mode = 

c^£^;s/c-o 7^X >i/i> J 
A - -A3(aK-077- . / n 

- — r A l- »/^ J 

Xv^c v» X 


Figure 11. Displacement on the antenna over a line 
through its center (origin). 
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Numerically, if Az = 1 m, D = 1010 m, X = 11.64 m, and X = 0.1224 m, 
then Ri 0-89 so tilt errors may be nearly as large as the 

subarray's beamwidth. 

We evaluate vibration in the x direction when 9^ = 0, 9^^ = 0, and 
when 

AZ = AZ cos(w^t) (37) 


where a)y/27r is the frequency of the vibrations. From (10), using (7), 
the signals e^"^^ and e^"^ in Figure 4(b) are 


e 


(+) 


e 


(-) 


^ ^ Sa[irX(ejj+8^^)/x] 

n=-N 

2ttX9 

• {cos[n( ■ y -- + 9q) + 


X 


+ cos[n(- 


27tX9, 


- 8q) - 






(38) 


For relatively small values of 9 the system's response will be 




However, the true response without vibration will very 
nearly equal the slope when 9^ = 0 times a„, so that the "error" in radians 

X A 

produced by vibration, denoted 9^, will be approximately 


= - [■ 


sub 


,(-) 




+ e J .= - [ 


2MX 




(39) 


2 M e' ' - ■ e 

after (27) is used. If normalized to the array basic beamwidth 9g we have 


9 

9 


V 

B 


,(-) 


2 e 


(+) 


while (38) becomes (for 9 = 0.5 9g) 

H ^ 


( 40 ) 
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A 

} = 0 


(-) 


N 

T ^ 

^ n=-N 


Sa[ 


M 0. 


(e +0 )] 

X xn'-* 


*{cos[n( 
+ cos[n( 


Ztt 0 
2ir 0 


X 


X 


M 

M 


+ 


+ 


A 

A 


(41) 


An examination of (41) shows that = 0 if 0^ = 0. Thus, for 
no antenna pointing error (0,, = 0) the presence of vibration produces no 

' ' A — — — » “ “* " 

error , regard! ess of the vi brati on mode , so 1 ong a_s jjt 2 ^ symmetri c . 

To study the error 0,,/Sp when 0 0, (40) was computed on a digital 

computer for various values of M, K, and 0 ^ when peak vibration dis- 
placement in each case is 1.0 m. Figure 12 illustrates curves of 0^ as 
a percentage of array basic beamwidth 0g versus the number M of subarrays 
in the array. Curves are shown for pointing error magnitudes |e /0 d| of 

X D 

0.05 and 0.1 for two vibration modes (K=l,2). The curves assume X=11.64 m, 
0=1012.63 m, and A=0.1224 m. It is clear that some combinationsof vibration 
mode and pointing error give rise to large errors due to vibration. 

For example, with M=7 subarrays in the linear array and mode 2 vib- 
ration (K=2), when pointing error 0 ^ is 5 % of the array’s basic beamwidth 
0g, the maximum error due to vibration of peak amplitude 1.0 m is 7.32% of 
0g. Figure 13 plots 0g versus M and gives 0g = 5.16min for M=7. Thus, 
maximum error due to vibration is 0.0732 (5.16) = 0.38min. Since the 
pointing error is 0.05 (5.16) = 0.268min, we find the vibration error even 
exceeds the pointing error to be corrected by the servo system. The vib- 
ration error cycles with time, and, for the example taken here, is illustrated 
in Figure 14. 
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Finally, we note that vibration will have an affect on signal -to-noi se 
ratio of the composite signals e^"^^ and e^”^. The affect on e^^^ is 
of great concern since signal -to-noise ratio for this signal should be 
maintained above 10 to 100 for good performance of practical normalizing 
circuits. Figure 15 shows the maximum loss in this signal-to-noise ratio 
for various M, K, and e^/0n when Tz ~ 1.0 m. 

X D 






3. RIGID BODY ANTENNA POINTING CONTROL 


In. this section, an antenna pointing control system is modeled and 
then designed assuming that both the solar collector and antenna are rigid 
bodies. The purpose of this design is two-fold. First, the design will 
generate technical data such as motor requirements, system response time, 
effect of noise and system parameter changes, and the stiffness of the 
control loop. Second, by examining time constant of the control system 
and the mode frequencies of the flexible collector, the effect of the 
collector vibration mode on antenna pointing accuracy can be estimated. 

The goal of the mechanical pointing control is to maintain a mechanical 
pointing accuracy of three arc-minutes. 

3.1 The Baseline Satellite Configuration 

Figure 16 depicts the baseline satellite configuration. The dimension 
of the solar energy collector and the position of the microwave antenna 
are shown in Figure 16(a). The size and the general shape of the yoke- 
antenna combination is shown in Figure 16(b). Figure 16(c) is included 
to show the size of the rotary joint connecting the yoke to the solar 
collector. 

The mechanical control of the antenna pointing is done by rotational 
motions about two axes. One rotational motion is about the Zy axis which 
is along the longitudinal center line of the collector as shown in Figure 
16(a). A second rotational motion is about the axis which is along 
the direction of hinges joining the antenna and the yoke as shown in 
Figure 16(b). 

Two coordinate frames are used in the sequel, one is the antenna 
frame (X^,Y^,Z^) and the other is the yoke frame (Xy,Yy,Zy). 
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(a) SPS satellite 



(c) The rotary joint (b) Antenna and yoke 


Figure 16. Baseline configuration for SPS satellite. 
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The physical parameters of the satellite are collected and listed 

here. For the collector: 

Mass: M = 16x10^ kg 
c 

Dimension = z - 17,000 m 
a = 4,200 m 
b = 600 m 

moments of inertia = 

(I.)? = ^ (a^+b^) = 19.5x10^^ kg-m^ 

c ^ 

(I L - = 1.54x10^^ kg-m^ 

C A ^ 

For the combination of antenna and yoke: 

Mass = = 13x10^ kg 

Diameter: D = 1,100 m 

Moments of inertia: 

(I^)z - |mR^ = " .983x10^^ kg-m^ 

(Ia^X^ - |-MR^ = .983x10^^ kg-m^ 


Other physical data include: 


slip-ring diameter = 350 m 

motor shaft diameter = .1 m 

♦ . . .. Antenna speed 

gean-train ratio a = — ^ — 

^ motor speed 


= 10 


-4 


The gear-train is used to connect the motor shaft to the rotary joint. 

the rotary joint has a diameter of 350 m and the motor diameter is .1 m 

-4 

a gear-train ratio of 10 is reasonable. 


Since 
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3.2 SHp-rinq Friction Characteri sties 

Conduction of electric power from the solar collector to the antenna 
yoke and from the yoke to the antenna is done by using slip-rings. The 
friction torque at each rotary joint is mostly due to the slip-ring friction 
torque which is a nonlinear function of speed of rotation. The nonlinear 
slip-ring friction characteri sties, as provided by NASA/MSFC, are shown 
in Figure 17. In the figure 0 represents the relative angular position 
between two sides of the rotary joint. The arrowheads indicate directions 
of motion. Notice that the friction torque has the saturation values of 
+ T^^„ = + 10^ N-M. In the region between the saturation values, the 
friction troque depends on e in a way resembling the compliance torque. 

c 

The compliance is the slope Kp which is -10 N-M per degree. The nonlinear 
friction torque exhibits the nature of shifting hysteresis. 

Figure 18 gives two analog electronic circuits which can be used to 
simulate the above friction characteristics on an analog computer. For 
system simulation using a digital computer, however, neither one of these 
circuits are convenient to use. To provide a convenient way for later 
digital simulation, a digital simulation flow-diagram is developed as 
shown in Figure 19. 

3.3 Estimation of Motor Parameters 

It is assumed that d-c motors are used as torques for antenna pointing 
control. The rating for the d-c motors are estimated here. Assume a maximum 
perturbation angle for the antenna to be [3] 

|a 0| =1 degree = .0175 radian 

' 'max ^ 

The moment of inertia of the yoke-antenna load is 
I^ = .983x10^^ kg-m^ 
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-4 

Given a gear-train ratio of a=10 , the reflected maximum angle perturbation 

and the moment of inertia of the load at the motor shaft is 


II = = 9830 kg-m‘ 

A A ^ 


A0 ' 



I 

a 


A0 


'max 


175 radians 


(see discussion in Section 3.6). 

The maximum slip-ring friction torque is 


= 10^ N-M at the ring 
= axlO^ = 100 N-M at the motor shaft 


The initial motor torque for a maximum step perturbation is 


' ^A I® Lax 
max 


U Imax . 1.72x10® 


N-M 


where T is the system time constant. Therefore the total maximum torque is 

.6 


T I “ Kr I 1 1 =100 + 

'max ' F'max ' 0‘max 


1.72x10^ 


(42) 


An estimate of motor power P is now given 


P = Average of [A0(t) x(t)] 


- Lax ‘"^Lax 


A0 ' I It I 
' max ' ' max 


1.75x10^ T^+3xl0® 


watts 


(43) 
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^ 23.46 T^+4.021x1Q^ ^p ^ 44 ^ 

Linder the nominal condition the steady-state motor power is to overcome 
the slip-ring friction torque. This power is 

Pg = HP (45) 

A list of P and t versus T is shown in Table 1. A baseline rating is 
chosen to be: 

P = 1681 watts = 2.253 HP 
T = 578 N-M 
T = 60 seconds 

Table 1. Motor Ratings 


T (sec) T (N-M) P (Watts) 


10 

17293 

301750 

30 

2010 

11694 

60 

578 

1681 

100 

272 

475 

300 

119 

69 

600 

105 

31 

1200 

101 

15 

1800 

100 

10 


Adopting the baseline motor rating, motor parameters are now estimated. 
Take the effective d-c voltage to be 


The average motor armature current is then 
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P 


1681 ^ 
120 


14 A 


" V 


dc 


The motor torque constant is 

.. - T . 573 . ^ N-M 

" I 14 ■ A 

sl 

Using the MKS system of units, the motor emf constant equals the motor 
torque constant in value, giving 

= 41.3 V-sec 

Assuming a 5% motor power loss, the motor armature resistance is computed 
by 

R = = .43f^ 

^ I- 


The result is summarized in Table 2 for the ease of future reference. 


Table 2. Motor Parameters 


Rated Power 

1681 

W 

Rated Voltage 

120 

V 

Torque constant, K.^. 

41.3 

N-M 

A 

Emf constant, 

41.3 

V-sec 

Armature resistance, R 

d 

.43 1 



3 . 4 Effect of Centrifugal Force Induced Torque 

Under the nominal condition the yoke-antenna assembly is rotating 
about the Zy axis of the SPS satellite at nearly constant speed. The 
rotation generates on each mass element of the antenna a centrifugal 
force which is directed perpendicularly to and away from the Zy axis. 
This force produces a torque about the axis of the antenna, which may 



39 

be viewed as a load disturbance torque for the pointing control system. 

This torque is analyzed with the help of Figure 20. 

Let AM be the incremental strip of antenna mass as shown in Figure 
20. Let be the nominal angular velocity of the yoke-antenna assembly 
about the Zy axis, which is very close to the Earth's rate of rotation. 

The incremental torque about the axis due to the centrifugal force on 
AM i s 

At^ = (AM)(acceleration)(torque arm) 

= (aX^ AZ^)(s^^ Z^ sin0)(Z^ cose) (46) 

where o is the area mass density of the antenna, which is assumed constant. 
Integrating over the entire antenna plane gives centirfugal torque as 

j-R 2 2 

''^c “ ^ ^ ^^A 

area 0 

2 T 0 

= 2 a sin2e J -Z^ Z^ dZ^ 

0 

where R is the radius of the antenna. After the integration the total 
centrifugal torque is obtained as 

T. = J a sin20 (47) 

Co c 

Maximum torque occurs when 0=45°. The Earth's rotation rate is 

-5 

= 15 deg/hour = 7.27x10 rad/sec 
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Therefore the maximum torque is 

l^cU = 2598.2 N-m (48) 

Comparing this torque to the maximum slip-ring friction torque, 
which is lO^N-m, the former is less than .3% of the latter. The effect 
of the centrifugal force induced torque therefore needs no special 
attention, but may be lumped into the slip-ring friction torque as a 
perturbation. It is expected that such a small perturbation will not 
cause problems for the antenna pointing control. Computer simulation 
will be used to verify this fact. 


3 . 5 Effect of Variation in Antenna's Moment of Inerita 

Since the antenna disc may rotate about the axis as indicated in 
Figure 21, its angular position e about this axis may change. As a result, 
the antenna's moment of inertia about the Zy axis may change. It is 
desired that such change will not have adverse effects on the accuracy 
of the pointing control. 

Let the moment of inertia matrix of the antenna expressed along the 
antenna frame (X^,Y^,Z^) be 


I 


0 

0 


i_ 


0 0 


Referring to Figure 21, the antenna frame and the yoke frame are related 
by an angle e. The antenna's moment of inertia expressed along the yoke 
frame (XY,Yy,Zy) is 
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where c9 = cos8 and se = sine. In particular, the moment of inertia about 
the Zy axi s is 

r = I S9^+I ce^ (50) 

z y z 

Note 

1 2 

I = o-M,R = max, ant. mom. of inertia 
y 2 A 

1 2 

I = X = niin* mom. of inertia 

z 4 A 

Taking the derivative of (50) with respect to 8 

~ = 2 se c9 - 2 se ce = (ly-I^)sin 2e 

Therefore the variation of the moment of inertia normalized with respect 
to 1 2 is 

- 1) sin 2e A8 = sin 29 a9 (51) 

For 9=45^, (51) is at its maximum. For a9=1 degree = .0175 radian, the 
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specified attitude control accuracy for the SPS satellite, 



“I 


—'max 


.0175 = 1.75% 


( 52 ) 


A plant parameter variation shown in (52) is too small to cause adverse 
effects in pointing control, and can be ignored in the design of pointing 
control. However, computer simulation will be performed to confirm this 
observation. 


3.6 Pointing Control System Modeling 

An all -attitude pointing control capability of the SPS antenna is made 

possible by allowing two degrees of rotational freedom for the antenna disc 

as shown in Figure 22. The first rotary joint allows the rotation of the 

antenna-yoke assembly with respect to the collector about the axis, while 

the second rotary joint gives rotational freedom for the antenna with 

respect to the yoke about the axis. It is assumed that the coupling of 

the antenna's motion about these two axis is negligible. This is due to 

the fact that the motion of the antenna is very slow and that the angular 

perturbati ons about the two axes are very small. Therefore the dynamics of 

the pointing control system for each rotary joint may be separately modeled. 

The First Rotary Joint Let 0» and e_ be the absolute angular positions 

c 

of the antenna-yoke assembly and the collector, respectively , about the Zy 
axis. Their dynamics are given by 

9a = T ' - Tf ( 53 ) 

I s^ 0 = -t'+Tj. (54) 

c c F 


I 
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where and are moments of inertia of the antenna-yoke assembly and 
the collector, respectively; x' is the control torque originated from the 
motor; and Xp is the frictional torque of the slip-ring. Rotational 
compliance and other kinds of friction are assumed negligible. Adding 
(53) and (54) gives 


Define the ratio of moments of inertia 

I, 


B = 




(56) 


Then (55) gives 


6'c = - 


(57) 


The relationship among the generated torque x, applied voltage V, and 
the shaft speed S0 for a separately excited d-c motor is 


X = K. 


V-KpS0 


T R 


(58) 


where K_, K- and R are defined in Section 3.3. Noting the presence of 
1 h a 

a gear-train with a ratio a, the following relationships exist. 


e = - ( 9 a- 0 ^) = - (1+6)0 

a A c a 


A 


(59) 


X = a X 


(60) 


Ignore the nonlinear friction torque Xp in (53) for the time being. Then 


'a «A = 


(61) 


Eliminating 6 ,x and x' from (53), (59), (60). and (61), yield 


Kt p p 

^V = [(a\)s2.-P(H6)s]^ 
a a 


( 62 ) 
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Let 


U = “ Ia 


(63) 


which represents the reflected moment of inertia of the antenna-yoke 
assembly seen at the motor shaft. The transfer function of the first 


rotary system including the motor is therefore 

a K, 

0, 


'A 




V 


K K 

s[5+t 1 D --(i+a)] 


'a 


(64) 


For later convenience, define 

®A " ct 


(65) 


which is the reflected absolute attitute of the antenna-yoke as seen by 
the motor. Using (64) in (65) gives 


®A 

Gp V 


Kj 


K K 

sls+jT— (1+e)] 


'A a 


( 66 ) 


Eq. (66) is the transfer function of the plant at the first rotary joint. 
Notice that the effect of the reaction motion of the finite mass collector 
is represented by B in the equation. For an infinite mass collector, B=0. 
A block diagram for this plant is given in Figure 23(a). 

In sections 3.1 and 3.3, values of various parameters have been 
obtained as: 

= 9830 kg-m^ = .43 Q. 

K-^ = 41.3 o = 10'^ 

Kg = 41.3 V-s 
















I 
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Quantity 6 as defined by (56) has a value of 


S 


. .983x10^^ 
IT 

19.5x10'^ 


.05 


With numerical values inserted, (66) becomes 

n - . .01 

" V " s(s+.424) 


(67) 


In this model the nonlinear friction was not included. The effect of 

the nonlinear friction will be treated as a torque disturbance in the control 

system design. 


The Second Rotary Joint For the control about the axis at the 
second rotary joint, the collector's moment of inertia may be assumed 
infinite. This is demonstrated as follows. Referring to Figure 24 and 
using the baseline data given in Section 3.1, moments of inertia about 
axis are: 

For the antenna: 

For the collector 

dc)x “ jMj. = ^(1.231M^)(31R)2 
A 

= 394.3 M 
= 1577 (1^),^ 


Thus, a rotation of antenna about X^ by 1° causes a rotation of collector 
about less than 3 arc -seconds. For a 3 arc-minute mechanical pointing 
accuracy, the assumption of infinite moment of inertia contributes only 
.11 arc-second error. 




51 

The following are further assumed: (1) the same type of d-c motor 
will be used at the second rotary joint as the type used at the first 
rotary joint. (2) Same gear-ratio. (3) Mass of yoke is small as compared 
to that of antenna. Under this condition the transfer function for pointing 
control at the second rotary joint can be obtained from (66) by setting s=0, 
yei Iding 


Q - ■ A 

S ■ T 


A a 


K K 

s[s+tt-4 ] 




( 68 ) 


where 9^ is the reflected antenna angle about the axis as seen by 
the motor. The actual antenna angle is 


♦a “ '•’A 


(69) 


With numerical values inserted, (68) gives 


A 01 

^p " T ~ s(s+.404) 

A block (diagram of the plant representing the seconci rotary axis 
in Figure 23(b). 


(70) 


is given 


3 . 7 Pointing Control System Design 
Design Background 

Under the assumption of negligible dynamical coupling between 
the two axes of rotation, the two associated control channels can be 
separately designed. The need of pointing control is to counter-act the 
effect of external disturbances and changes in plant parameters. External 
disturbances include gravity gradient torque, solar pressure torque, and 
collector motion coupled through the rotary joint. Change in plant para- 
meters include the change in the antenna's moment of inertia due to change 
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in antenna attitude, change in motor armature ressistance caused by themial 
conditions, aging effects, change of structural shapes and dimensions 
due to temperature changes, and others. Most of the disturbances are corr- 
stant or very low frequency cyclic. Parameter changes are also very low 
frequency in nature. It is assumed that the frequency range for distur- 
bances and parameter changes is below .001 hertz (15 minutes per cycle). 
Effects of instrument noise will be ignored because they will be well 
attenuated by the low-pass characteristics of the system. 

Linder the nominal condition, the antenna-yoke is rotating at the 
first rotary joint with a constant rate of one. revolution per 24 hours. 
Magnitudes of disturbances at this joint are assumed small so the direction 
of rotation is not reversed. The second rotary joint is stationary under 
the nominal condition. With disturbances added, it may rotate in either 
di recti on. 

Comparing (67) and (70) shows that they differ only in time constant 
of an amount of less than 5%. Many well designed control loop for one plant 
should therefore perfonn well for the other. In view of this, a control 
system design will be made for (.67), the plant for the first rotary joint. 

Figure 25 depicts an antenna pointing control loop for the first rotary 
joint. An ideal gear-train is assumed. The pointing error is sensed 
by the microwave error sensing method presented in Section 2.1. An 
algorithm is used to generate, from the error signal, a command signal 
for pointing control. Three ways of generating command signal have been 
developed previously. [1] For the convenience of design the transfer 
function between error e and command signal may be taken to be unity. 

The present control system design is to determine the block labeled "network 
and amplifier gain". 



r 


Network 

and 

ampl ifier 
gain 


u 


Command 

signal 

generation 


Figure 25. 


Antenna 




Gear-train 


the pointing error 


pointing control for the first rotary axis. 


CO 
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The Design Except for the 3 arc -minute mechanical pointing accuracy 
requirement, no other system performance specifications were given. 

Linear quadratic optimal control design will be adopted. Since the plant, 
represented by (67), has single-input and single output, frequency domain 
optimization by way of the root-square-locus technique v/ill be used. [4], 
Because of the constant rate rotation under the nominal conditions, the 
control system for the first rotary joint should be a follow-up system. 
Follow-up design will also be used for the second rotary joint to provide 
flexibility for the possible need of follow-up operations. 

Begin with (67), the plant transfer function* 

r - -01 /In 

■ s(s+.424) 

which has poles at p-j=0, P 2 = 0 » and P2=--424. Open-loop poles for the 
root-square locus plot are 



There are no open-loop finite zeros. The root-square loci are shown in 
Figure 26. Choose the closed-loop root-square poles at 

& Q 2 = .8 + j 1.48 
Q 3 = -.iS-q^-Q^ = - 1.78 

Then 

-Q^ & -Q^ = -.3 + j 1.48 
-' 53 = 1-78 

ic 

An integrator is added to make the system type-2. 
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The closed loop poles for the control system are given by 


q = - /-Q ^ 1.297 1 120.8 


Compute 


qo = - AQo = 1.297 1 -120.8 




q^q2 == 1 . 682 
q^j+q^ ‘ -1.328 

= - .79 

q^Pa 


The reciprocal of the closed loop zero is 


T = . (i . X . X) = 

‘’l- ^2 ^3 


PiPa pj' 


= .79 + .75 = 1.54 


The closed- loop transfer function is 


H = 


- ®A . (1+Ts) 


IT 




l+Ts 


,, s,,, '’r'^2 , 1 2, 

(1 )(1 s + s ) 

q/ . q^q9 ^^^2 


1 + 1.54 s 


(l + .75s)(l + .79s + .5955"^) 


1 + 1.54s 


1 + 1.54s + 1.18s^ + .45 s^ 


( 71 ) 


( 72 ) 
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Step response simulation was performed to test (72). The overshoot 
was found to be more than 40%: not a good response. Repeated trial- 
and-error adjustment led to the following closed-loop transfer function. 




1 + 3.164s 


1 + 3.164s + 3.082s^ + .595s^ 


This system had a step responce overshoot of 18.7% and a time-to-peak 
of 3.4 seconds. The corresponding open-loop transfer function Gg was 


obtained from 


giving 


1= 1 + i- 

H Gg 


G = .3245(1 + 3.164s) 
° s^(l + .193s) 


(73) 


(74) 


(75) 


Two compensation configurations were used as shown in Figure 27; 
the first is a series compensation and the second is a series compensation 
with rate feedback. The second configuration offers better counter-action 
for load disturbance and plant parameter variation. 

For series compensation: 

„ . (1 + 3.164s)(1 + 2.36s) , 

s(1 + .193s) ^ 

For series compensation with rate feedback: 


choose K.J. = 475.76 - 476 


then 6^ 
c 


168,1(1 + 3.164s) 
s 


(77) 

(78) 


A block diagram for the simulation of the designed closed-loop 
system is given in Figure 28. The dashed line part is used only for 
series compensation with rate feedback. Figure 29 shows the step response 
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Figure 29. Step response of the nominal system 
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of the system under the nominal conditions. The reference input is 175 
radians for motor, corresponding to 1° for the antenna. It should be 
mentioned that the system is not expected to be subjected to a step 
function type of corrmand. Step response simulation is used merely to 
reveal the system's transient characteristics. 

System Sensitivity The response sensitivity of the designed system 
was studied by way of computer simulation. Results are given in Table 
3. As a whole, the response did not change much for parameter change 
as high as 20%. Also, the superiority of the series compensation with 
rate feedback is evident. 

The Constant Drag Torque Since the first rotary joint does not change 
its direction of rotation the nonlinear slip-ring friction torque appears 
as a constant drag torque with a magnitude of entering the system 

at the point as shown in Figure 30. Results of a step response simulation 
are given in Table 4. One sees that the control is capable of nulling 
the effect of constant drag torque very effectively. A CSMP listing of the 
simulation is attached as Appendix A. 


Table 4. Response of System with Prig Friciton 


Overshoot n 

Time to peak tp 

Error at 4 t : 

P 


motor angle 
antenna angle 


21 % 

3.3 sec 
.13 rad 
.04 arc-mi n 


Theoretically, the steady-state pointing error due to a constant 
drag torque should be zero. From Figure 30, the response of the system 

fl 

output due to a constant disturbance torque 



r 

System 

condition 

Series 

compensation 

1 

Nominal 

K,/R 10 % up 

1 d 

10 % down 

" 20 % up 


" 20 % down 

Series 

compensation 
with rate 
feedback 

Nominal 

” % up 

" 10 % down 

" 20 % up 

" 20 % down 


Overshoot 

n 

T ime 
to peak 
t 

P 

Settl i ng 
time 
t 

• s 

19 % 

3.4 s 

9.0 s 

16 % 

3.2 s 

9.6 s 

21 % 

3.6 s 

8.8 s 

15 % 

3.2 s 

9.8 s 

24 % 

3.8 s 

8.6 s 

■ 

19 % 

3.2 s 

8.8 s 

19 % 

3.2 s 

8.8 s 

20 % 

3.2 s 

8.8 s 

18 % 

3.2 s 

8.8 s 

20 % 

3.2 s 

8.6 s 











































64 


9 ^( 3 ) = 


1 


( 79 ) 


+ ^CG(.(s) + Kg(l + 8)s] 

3 


The steady-state value of 0^ is 


8^ (t->co) = s0^(s^O) = 0 (80) 

where the property G^(s->0) = °° has been used. 

The Nonlinear Friction Torque The second rotary joint has little 
rotation under the nominal conditions. To counter anamalous effects, 
it may rotate in either direction. Thus the shifting hysteresis type of 
nonlinear friction torque, as portrayed in Figure 17, may come into play. 

This friction torque is a nonlinear function of the speed of the motor, 
and it enters the system in a way as shown in Figure 31. The nonlinear 
function is designated Gp. Gp can be simulated according to the flow 
diagram of Figure 19. Note that all load quantities have been reflected 
into quantities as seen by the motor; the value for Kp in Figure 19 should 
be handled similarly. Using the data in Figure 17, 




- 10 


6 ^ 
Deg 


at antenna 


57.3x10 

10 ^ 

l7^ 


N-M 


rad 


= - .573 


N-M 

rad 


at motor 


max 


10^ N-M 
10 ^ 


at antenna 


10 


100 N-M 


at motor 



Centrifugal force 
induced torque 

'c- 


350 rad. 


Nonlinear friction torque 



nonlinear friction torque 
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To reveal the effect of nonlinear friction torque, a square-wave may be 
used as reference input for system simulation since rate reversal is 
needed to exhibit the shi-fting hysteresis characteristic of the nonlinearity. 
In order to drive the nonlinear friction torque to its saturation level 
the amplitude of the square-wave should be greater than 1° in terms of 
the actual antenna angle. An amplitude of 2 ° is chosen. The corresponding 
amplitude for the motor angle is (2x10 )/57.3 = 349.04 radians. A simulation 
was performed in the manner described. Results of simulation are listed in 
Table 5, which shows the adequancy of the control. Appendix B contains 
the CSMP listing of the simulation. 


Table 5. Square-Wave Response of the System having 

Shifting Hysteresis Type of Friction Torque 


Overshoot n 

Time to peak tp 

Error at 4 tpi Motor angle 

Antenna angle 


20 % 

3.19 sec 
. 1 4 rad 
.05 arc -min 


The Centrifugal Force Induced Torque Recall from the discussion 

in Section 3.4 that the maximum torque on antenna caused by the centrifugal 

force is less than .3% of x which is 100 N-M as seen by the motor. Insert 

max 

T = 1 N-M, which is 1% of i into the last simulation program and repeat 
c m3x 

the simulation. The results are nearly the same as those in Table 5, 
showing that the pointing control is capable of nulling the effect of x^ 


satisfactori ly. 



3.8 Effect of Noise in Sensed Pointing Error Signal 

It has been shown that, based on the set of given assumptions, the 
designed mechanical pointing control system more than adequately meets 
the 3-arc-minute requirement for pointing accuracy. If the microwave 
sensed pointing error is noise contaminated, however, the situation 
could be quite different. As discussed in Section 2, the noise in sensed 
error is low frequency in nature, which can be as low as d-c. Under this 
condition the noise directly dictates the accuracy of the overall pointing 
control system. At this moment, it is not practical to arrive at a 
numerical value for the noise level. The remark here serves to emphasize 
the possible adverse effect of the noise. 

3 . 9 E ffect of Collector Vibrations 

It has been discussed in Section 3.7 that the frequencies of vibration 
for the collector is less than .001 hertz which is considerably lower than 
.2 hertz, the approximate bandwidth of the pointing control system. Thereto 
one expects errors caused by the collector vibration be effectively 
diminished by the pointing control. To support this claim, computer sim- 
ulations were performed as follows. 

Pointing error caused by collector vibration is equivalent to the 
error of the closed-loop control system excited by a sinusoidal reference 
input. Reference inputs having frequencies of .001 hertz (about 16.6 
minute/cycle) and .01 hertz (about 1.66 minutes/cycle) were used in the 
simulations. Results are shown in Table 6, which confirm the claim. 
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Table 6. Tracking Error Due to Collector Vibrations 


First simulation 


vibration frequency = 
vibration amplitude = 
maximum steady-state error = 


.001 hertz (=; 16.6 min/cycle) 
1 degree 

.0055 rad. at motor 
,002 arc-mi n at antenna 


Second simulation 

vibration frequency = 
vibration amplitude = 
maximum steady-state error = 


.01 hertz {- 1 . 66 min/cycle) 
1 degree 

. 55 rad , at motor 


.19 arc -min at antenna 



4. APPROACHES FOR MODELING THE FLEXIBLE COLLECTOR BODY 


It is recognized that the structural dynamics of the flexible collector 
body may be modeled by transforming it to an equivalent orthotropic plate. 
Then the wealth of well known theory of plates can be applied readily. 
However, a thorough study of the structural responses of the collector 
plate requires a major effort. The discussion here will focus on the 
two approaches whrich may be used to obtain the desired model. They are the 


analytical approach and the simulation approach. 

The Anal y tical Approach In this approach the structure of the collector 
truss is first transformed to an equivalent grid, which is then trans- 
formed to an equivalent orthotropic plate as shown in Figure 32. The 
differential equation for an orthotropic plate, including the inertial - 
force term, has the form [5] 


3 ^ 


0 + 2H 


3x^9y^ -^sy^ at 


(81) 


where w = u)(x,y,t) is the deflection in the z-direction as shown in 

Figure 32(c), q is the external loading per unit area, h is the thickness, 

and p is the density of the plate material. The problem of modeling, simply 

stated, is the determination of constants D , D , and H. 

x’ y* 

Transformation from grid work to an orthotropic plate has been 
thoroughly studied by Timoshenko and his cohort [1]. Let the geometry 
of a grid be def lined in Figure 33, then the constants for an equivalent 
orthotropic plate are given by 


D.. = 


E h' 
12 


D = 


E h' 
12 
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(a) The collector truss 





(c) An equivalent orthotropic plate 


Figure 32. Transformation to an equivalent orthotropic plate 
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( 82 ) 


where E is the modulus of elasticity, and c-j and C 2 are torsional 
rigidities of the beams parallel to the x and y axes, respectively. When 
the beams have rectangular cross-section as shown, c^ and c^ are given by [6] 


_ G ,3 . n 192 "l y 
3 ^1 " 5 ^ 


C2 = I h[l 


4r tanh ~) 

n=l,3,5— 


(S3) 


192 2 7 1 . u ,nTT h 

~ li" n-l ? ^ ‘2' ^ 

IT n- i , 3 , n c 

5.— 


For most practical cases the infinite series converges rapidly; usually 
two or three terms are sufficient to give accurate result. For example, if 

4^ = 1 , then 
^1 

I -4-ir tanh ^ 

n=l,3,5— h^ ^ 


= tanh 5- + 4=- tanh 4^ + tanh %^ + - - - 

2 36 2 56 2 

= 0.91715 + ^ (0.9998) + ^(D + 


= 0.91715(1 + .004486 + .0003489 + ) 


= 0.9216 


Similarly, the results of the series become 1.00073 and 1.00414 for 

h h 

equals 2 and 4 respecti vely . Using the result for — = 2 in (83) gives 

ti t-, 



r 


73 

c-| ~ y (y) -00073)] 

IT 

= I (-fp)^ (0.6861) 

= I (j) (0.6861) 

= 0.0286 G h'^ . 


Clearly C 2 =c-| if bending moments and twisting moments are 

related to the deflection 03 by the following expressions. 


t, .2 

Eh 1 9 0) 

X ■ " 12 b. 2 

I 9x 

Eh^ h 9^(j 

y == “ ITiTTT 



I a 03 


1 9x9y 


(84) 
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Analytic methods for transforming a truss structure as shown in Figure 
32(a) to an equivalent grid, or directly to an equivalent orthotropic 
plate have not been developed. It is this part of the study which requires 
a major effort. However, some related cases have been studied, such as 
plates with many holes [7-10] and corrugated pipes. [7-11]. The 
equivalent solid plate or smooth pipe will produce the same strain, defor- 
mation, or strain energy as the actual perforated plate or corrugated 
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pipe. The results for the actual structure have been carried out by testing, 
numerical calculation through point-matching or collocation technique, 
and other appropriate analysis such as energy method. 

The Simulation Approach Constant D^, D^, and H of Eq. (30) for a 
given truss may be obtained by a pa^^ameter matching technique with the 
help of computer simulation. Well known computer programs exist, such 
as NASTRAN and STROYN, which use finite element approximations for 
flexible structure. 



5. CONCLUSION 


A study of the pointing control system for SPS was accomplished. The 
study included two major areas, namely, the analysis of a pointing error 
sensing method for the SPS antenna, and the design of a mechanical pointing 
control system for the antenna. In addition, approaches for modeling the 
flexible body collector dynamics were given. 

In the area of point error sensing, an analysis was made on a pre- 
viously proposed method for developing pointing angle signals to be used 
by a servomechanism to point the antenna toward a rectenna on ground. From 
the analysis result, the method was found to be feasible. A normalized 
technique was added to the method to produce error signals that are in- 
dependent of the received signal magnitude. Effects of vibration and noise 
on pointing control accuracy were analyzed and presented in the form of 
charts . 

In the area of mechanical pointing control of the antenna, a baseline 
satellite configuration was selected as shown in Figure 16 (p. 30) together 
with a set of baseline parameters. Analysis were made on various phenomena 
which might influence the accuracy of the pointing control. They include 
the nonlinear slip-ring friction of rotary joints, motor characteris tics , 
effect of centrifugal force induced torque, effect of variation in antenna's 
moment of inertia, effect of finite collector mass on pointing control, and 
the effect of vibration modes of the collector. The friction and backlash 
character! sties of the gear train was not available for analysis and were 
ignored. A mathematical model of the plant, whose output is the antenna 
pointing direction, was developed. The pointing control system was designed 
using an optimal control approach with the help of the root-square-1 ocus method. 
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The mechanical pointing accuracy, with the various anomalies present, was found 
to be well within 1 arc -minute. 

From the results of this study it appears that the desired pointing 
accuracy of + 3 min is possible. Suppose a reasonable error budget is 
+ 1 min for each of the major error sources of (1) noise, (2) vibration, 
and (3) mechanical effects (gear backlash, stricktion, etc.). Suppose 
further that we require the + 3 min to correspond to at most 50% of the 
beamwidth 6g so that operation is always in the linear part of the error 
pattern. Then from Figure 13 and 0g = 12 min we require M<3 subarrays in 
the linear array. The allowed vibration error is therefore 1 min/ 12 min - 
0.083 or 8.3% of 9g maximum. From Figure 12 this allowance is not exceeded 
for M=3 even for mode 2 (k=2) vibration (az = 1 m). The noise error of 
+ 1 min can be achieved by providing sufficient power in the ground station 
pilot transmitter. 
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APPENDIX A. CSMP FOR THE FIRST ROTARY JOINT CONTROL SYSTEM. 


«:(c:{::}:CQNriNUaUS SYSTEM MODELING PRGGRAM=J'*=J'* 

«=i^«PRQBLEM INPUT STATEMENTS’!'^^ 

STEP RESPONSE WITH DRAG DISTURBANCE 
R=175, 

TD= I 00,0 
E=R-C 

X1=13.754'E 
X2=INTGRL( 0.0, XI ) 

X3=LEDLAG( 3. 164, 0. 193 ,X2) 

X4=DERIVI 0.0,X3J 

X9=2.36*X4 

V=X9+X3 

E3^V-C2 

C2=43.4’i'Cl 

C1=OERIV(0.0,C) 

C2 = 43.4=i‘Cl 
X5 = 96, OS^f'ES 
X6=TD+X5 

X7=0. 00010173«X6 
X8^INTGRL(0-0,X7) 

C=INTGRL (0.0, X8) 


TIMER FINT1M=20., GUTDEL=.2, DELT=.001, PR0EL=.2 
PRTPLT £, V, C 

LABEL STEP RESPONSE WITH DRAG DISTURBANCE 

END 

STOP 
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APPENDIX B. CSMP FOR THE SECOND ROTARY JOINT CONTROL SYSTEM. 


t.*^:*:CQNTINUQUS SYSTEM MODELING PROGR/iM=s^={=^« 
***PRCBLEM INPUT ST AT EMENT 

« POINTING CONTROL WITH SHIFTING HYSTERESIS TYPE 

* OF NONLINEAR FRICTION TORQUE 

# KF=--573, D£LT=,001» D=KF*D£LT 
INITIAL 

CONSTANT SGNP=0. rTM= 1 00 . t T DP=0. ,TD=0, ,D=-. 000573 
DYNAMIC 
SORT 

R = 3 5 0-«STEP( 0,0)-200.=«'STEP (10,0) 

E = R-C 

Xl=25, 18«E 
X2=INTGRLI O.OtXI) 

X3 = L ED L AG ( 2,20 ,0.270, X2 ) 

X4-DERIV (0.0, X3) 

X9=2. 36*X4 
V=X9+X3 
E3^V-C2 
C2=43.4^C1 
C1=DERIV(0-0,C ) 

C2=43.4=^^C1 
X5 = 96. 05>i=E3 
X6=T0+X5 

X7=0. 000 10173>4=X6 
RATE= INTGRL(0.0,X7) 

C= INTGRL ( 0. 0, RATE) 

NO SORT 

SGN=S IGN ( 1. 0, RATE) 

IF (SGN .EQ. SGNP) GOTO 301 
TDP=0. 

302 TD=TDP+0*RATE 

IF (A3S(TD) ,GT, TM) TO=S I GN ( T M, TD ) 

GO TO 304 

301 IF (ABS(TDP) .EQ. TM) GOTO 303 
GOTO 302 

303 TD=TDP 

304 TDP=TD 
SGNP=SGN 

SORT 

TERMINAL 

TIMER FIMTIM=20., QUTDEL=-2, DELT=.001, PRDEL=.2 
METHOD RKSFX 
PRTPLT C, E, RATE, TD 
LABEL STEP RESPONSE 
END 
STOP 
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